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I.  INTRODUCTION 


A  previous  report^  proposed  a  drag  determination  procedure  for 
hypervelocity  "Kinetic  Energy"  projectiles  intended  for  anti-tank  application. 

A  similar  estimating  procedure  for  use  in  antiaircraft  fire  for  spin 
stabilized  projectiles  is  now  offered.  The  basic  boundary  conditions  are 
retained,  i.e.,  sea  level  air  density  and  flat  trajectory,  and  it  is  assumed 
that  the  projectile  geometry  is  reasonably  conventional. 

The  calculation  is  quite  similar  to  that  of  Reference  1,  except  that  the 
fin  contribution  is  not  required.  Also,  retardation  and  time-of-f light  print 
outs  have  now  been  included.  The  results  are  presented  in  desk-calculator 
form  and  the  program  is  included  as  an  Appendix. 

II.  PROCEDURE 

Figure  1  is  a  schematic  of  a  typical  spin  stabilized  projectile  and 
Figure  2  defines  the  nomenclature  employed  in  the  analysis.  Standard  sea 
level  air  properties  are  assumed. 

A.  A  conventional  partition  of  the  drag  coefficient  as  composed  of  wave, 
viscous  and  base  components  is  assumed.  The  wave  drag  coefficient  is  taken  as 


CDW  * 7  M 


-.28  -  -1.73 


directly  from  Reference  2, which  examined  a  large  number  of  spin  stabilized 
projectiles  to  the  Mach  5  region  and  correlated  the  experimental  data  to 
develop  estimating  criteria.  Figure  3  indicates  the  transposition. 

The  hypersonic  base  drag  coefficient  is  found  by  inference  from 
experimental  data  on  cones.  A  patch  procedure,  as  described  in  Appendix  A,  is 
imposed  and  the  result  is  a  bilinear  characteristic  from  M=2  to  M=10  with  the 
knee  at  M=5.  Available  references  offer  little  insight  into  the  aerodynamics 
of  hypersonic  flow  in  the  wake  of  cylindrical  bodies.  The  wake  flow  behind 
cones  has  been  investigated,  however,  and  the  results  of  these  open  literature 
studies  is  included  in  Appendix  A.  Thus, 

CDB  =  (.050  -  .0034M)  db2 


is  used  for  use  in  the  range  5<M<10. 


W.F .  Donovan ,  "Hypothetical  Zero  Yaw  Drag  Coefficient  of  Kinetic  Energy 
Projectiles  Between  M=5  and  M=10 ARBRL-MR-03041 ,  August  1980,  ADA  it-  090009. 

W.F.  Donovan  and  Susan  A.  Wood,  "Automatic  Plotting  Routines  For 
Estimating  Properties  Of  Spin  Stabilized  Projectiles  In  Flat  Fire 
Trajectories  At  2<M<5,"  ARBRL-MR-03204 ,  October  1982 ,  AD  #120658. 
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Figure  1.  Typical  Spin  Stabilized  Projectile 
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Figure  2.  Projectile  Nomenclature 
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Figure  3.  Nose  Wave  Drag  Coefficient  Correlation 
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The  viscous  drag  component,  C^y,  is  obtained  by  Mach  extrapolation,  Figure 
4,  of  the  friction  factor  as  employed  in  Reference  1  and  of  the  effect  of  the 
ogive  as  empirically  found  in  Reference  2. 


Cm7  =  • 000173  (13.84  -  1.184M)  (.083  A  +  .0625  i  )  (30)  EXP 

an 


(4.6/R) 

3 


used  for  the  range  5<M<10.  The  viscous  flow  mechanics  have  been  examined  from 
divergent  assumptions  (Refer  to  Appendix  B)  and  are  considered  here  from  the 
most  conservative  viewpoint. 

For  a  projectile  with  a  supercaliber  rotating  band,  Reference  3  suggests 
a  constant  .015  C^R  in  the  lower  velocity  regions.  This  value  can  be  expected 
to  decrease  in  the  hypervelocity  regime  and  is  taken  here  as 


CDR  =  . 07143M  dr. 


The  total  drag  coefficient  is  then  the  sum  of  its  parts,  or 


CD  CDW  +  CDB  +  CDV  +  CDR* 

B.  The  Mach  number  along  the  trajectory  is  given^  as 


M 


where 


M  =  Mach  number  at  range  ns,V 
J  =  Operational  Parameter  =  k 

Q  =  Operational  Parameter  =  p 


2m  * 


s  =  Range, 

b  =  Intercept  of  vs  M  characteristic, 
k  =  Slope, 

Aref  =  Reference  area,  and 
m  -  Mass  of  projectile. 


L.Mm  Freeman  and  R.H .  Korkegi ^  " Projectile  Aft-Body  Drag  Reduction  by 
Combined  Boat-Tailing  and  Base  Blowing fr  AF  API-TR-7  5-111 3  February  1976 . 

WJ?.  Donovan >  nSimplified  Determination  of  Retardation  For  Kinetic  Energy 
Projectiles /'  BRL  Memorandum  Report  No.  ARBRL-MR-03020 ,  May  1980  y  AD 
# 086095 . 
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Figure  4.  Conversion  of  Mach  Number  to 


C.  The  average  velocity  decrement  over  the  selected  range  is  determined 
by 


AV 


(v 

v  sonic 


). 


where  AV  =  the  average  velocity  decrement  over  range  "s"  and  vson£c  =  ambient 
sonic  velocity. 


D.  The  time  of  flight  is  determined  from  the  average  velocity  during  the 
flight . 


TOF 


_ 2 _ s _ 

(M  +  M)  (V  .  ) 
0  sonic 


E.  The  input  parameters  required  by  the  program  are  listed  in  Table  1. 
Entries  to  the  HP-97  calculator  storage  are  in  caliber  units  where  the 
reference  diameter  is  either  inches  or  millimeters.  The  selected  range  is 
given  as  meters.  The  projectile  weight  may  be  inserted  as  pounds,  grams  or 
cubic  calibers.  It  is  necessary  to  employ  the  secondary  storage  for  recurring 
constant  numbers.  They  are  listed  in  Appendix  C. 

TABLE  1.  INPUT  PARAMETERS  FOR  PROGRAM  OPERATION 


HP-97 


REGISTER 

SYMBOL 

IDENTIFICATION 

UNITS 

A 

B 

C 

D 

dr 

Rotating  Band  Diameter 

Cal 

E 

s 

Range 

Meters 

I 

5 

Lower  Mach  Limit 

0 

M0 

Muzzle  Mach  Number 

1 

SL 

■Q 

Length  of  Nose 

Cal 

2 

Length  of  Afterbody 

Cal 

3 

db 

Base  Diameter 

Cal 

4 

R 

Ogive  Radius 

Cal 

5 

-  m 

Mass  of  Projectile 

Gram 

6 

m 

Mass  of  Projectile 

Pounds 

7 

m 

Mass  of  Projectile 

Cal3 

8 

d 

Diameter  of  Projectile 

mm 
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III.  RESULTS  AND  CONCLUSIONS 


For  the  projectile  described  in  Figure  1,  the  input  parameters  for  Table 
2  are  determined  and  the  calculator  program  prints  out  the  results  as  shown  in 
Table  3. 

Since  no  range  data  is  available  for  comparison,  verification  is 
lacking.  However,  within  the  regime  where  cone  results  are  applicable, 
indications  are  that  the  procedure  is  reliable  and  usable  for  preliminary 
estimating.  Where  alternate  formulations  of  the  drag  coefficients  of  the 
components  (wave,  viscous  and  base)  can  be  established,  their  substitution 
into  the  program  is  quite  simple.  The  subsequent  operations  are  identical. 


TABLE  2.  INPUT  PARAMETERS  FOR  HS-831  PROJECTILE 


ENTRY 

REGISTERS 

SYMBOL 

VALUE 

UNITS 


A 

B 

C 

D 

E 

1 
0 
1 

2 

3 

4 

5 

6 

7 

8 

9 


M  (initial) 


Mr 


n 


R 

m 


1  .06 
1500 
5 
7 

2.81 
1  .90 
.94 


12  .57 
361  .35 


d 


29  .9 


Cal 

Meters 


Cal 

Cal 

Cal 

Cal 

Grams 


mm 
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TABLE  3.  SAMPLE  PROGRAM  OPERATION 


Primary  Registers 


Secondary  Registers 


7. 666688666 

6 

25. 46888886 

O 

2. 8 i 6888888 

i 

8.686284 886 

4 

1 . 968888866 

z 

35.37888886 

£ 

6. 948668886 

8.688817344 

3 

12.57888886 

*T 

8. 888888888 

4 

361.3588866 

£- 

8. 682392886 

0 

8. 6888 88686 

0 

8.836111686 

£ 

6. 668686866 

l 

- 8.883488888 

»’ 

25.98886886 

O 

6. 68 3888888 

S 

8. 668888666 

r 

6. 668888888 

> 

8. 688888866 

ti 

8. 668886886 

B 

8. 668888888 

w 

1.668888886 

B 

1568.888866 

E 

5. 888888866 

i 
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Program  Output 


C 

M 


B 


Range  in  calibers _ _ _ 

Mach  number  at  range  _ _ _ _ 

Retardation  in  Meters/sec/kilometer 
Time  of  flight  in  sec  _ _ _ _ 


0.  1  26662168 

44  + 

5. @00000000 

4+  + 

8.118723975 

4+4 

6. 000000000 

44  + 

0. 182976366 

4+4 

7. 000000000 

44+ 

0. 896182885 

4++ 

8. 0000000t?0 

444 

0.09006004] 

44  + 

9. 000000000 

4++ 

0. 8644299o£. 

**4 

1 6.00000066 

++  + 

58167. 12376 

+  44 

6.861768514 

44  + 

213.9186826 

+++ 

0. 652494549 

4+  + 
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APPENDIX  A 


PATCH  PROCEDURE  FOR  TRANSITION  TO  HYPERSONIC  REGIME 

The  base  drag  coefficient  for  a  square  base  projectile  in  the  range  2<M<5  is 
given  by 

CDB  =  (.133  -  .02  M) 


and  acquires  a  value  of  .033  at  M-5.  For  an  assumed  decrement  to  .016  at 
M=10, 


.133  -  .02  M  =  -  k  M  +  b 


or 


.033  =  b  -  5  k 

with  .016  =  b  -  10  k. 


k  =  .0034,  and 

b  =  .05. 


Thus  , 


CDB  =  .05  -  .0034  M. 


The  skin  friction  coefficient  is  similarly  determined.  At  M  =  10  the 
extrapolated  decrement  produces  Cp,  =  .0002.  This  leads  to 


28.75  -  4.166  M  =  z  -  w  M, 


where  z  and  w  complete  the  linearization  in  M. 


Then  , 


z  =  13.84 


21 


with 


w  =  1 .184 

are  suitable  coefficients  for  the  range  5<M<10 . 
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APPENDIX  B 


DISCUSSION  OF  BASE  DRAG  COEFFICIENT 


T}_  1  TJ _ O 

Lyons  and  Brown  and  Zarin  offer  results  of  work  on  cones.  The  Lyons  and 
Brown  base  drag  coefficient  assumes  a  perfect  vacuum  in  the  immediate  wake  of 
the  body  while  the  Zarin  data  are  predicated  on  pressure  measurements  on  the 
model  mounted  in  a  wind  tunnel  facility.  McCoy^"~^  calculates  base  drag  at 
lower  Mach  numbers  for  flight  bodies  on  the  basis  of  a  Prandtl-Meyer  expansion 
around  a  sharp  corner.  Extrapolation  to  the  higher  Mach  numbers  precisely 
duplicates  the  Lyons  and  Brown  data.  The  results  are  compared  with  values 
from  this  report  in  Figure  B-l . 


Figure  R-l  . 

The  linear  assumption  of  the  present  report  intercepts  the  Zarin  data  at  M  =  6 
and  coincides  with  the  Lyons  and  Brown  -  McCoy  analyses  from  M  =  8  to  M  = 

10.  The  viscous  contribution  to  the  drag  is  treated  by  Lyons  and 


W.C .  Lyons 3  Jr.  and  H.S.  Brown 3  " The  Drag  of  Slightly  Blunted  Slender 
Cones  "  NOLTR  68-3  3  January  1968. 

N.A .  Zarin 3  " Base  Pressure  Measurements  On  Sharp  and  Blunt  9  Degree 
Cones  At  Mach  Numbers  From  3.50  to  9 .20  f1  BRL  MR  1709 3  November  1965 3 
AD#  369084 . 

Robert  L.  McCoy  3  ”MC  Drag  -  A  Computer  Program  For  Estimating  The  Drag 
Coefficients  of  Projectiles f*  ARBRL-TR-02293 3  February  1981 3  AD# 
A098110. 
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Brown  as  a  boundary  layer  phenomena  with  additional  components  due  to  induced 
pressure  and  transverse  curvature  effects.  Zarin  considers  the  viscosity  to 
be  negligible  in  comparison  with  other  terras.  Stivers  offers  a 
conventional  treatment  whereby  the  laminar  regime  is  superseded  by  a 
transitional  and  turbulent  flow,  and  then  converts  the  body  of  revolution  to 
equivalent  flat  plate  configuration.  This  present  report  simply  extrapolates 
from  lower  Mach  number  data  to  estimate  the  high-end  base  drag  coefficient. 
The  results  show  qualified  agreement  with  the  Lyons  and  Brown  argument. 


L.S.  Stivers,  Jr.,  "Calculated  Pressure  Distributions  and  Components  of 
Total  Drag  Coefficients  For  18  Constant  Volume  Slender  Bodies  of 
Revolution  At  Zero  Incidence  For  Mach  numbers  From  2.0  To  12.0  With 
experimental  Aerodynamic  Characteristics  For  Three  Of  The  Bodies ,"  NASA 
TN  D-6536 ,  October  1971. 
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LIST  OF  SYMBOLS 


Aef 


"DB 


"DR 


"DV 


"DW 


CF* 
Cpt  T 


D 

J 


M 


Q 


R 

Re 

TOF 

b 

d 


Area  of  cross  section  of  projectiles 

o 

Reference  area  (  .785  cal^) 

Total  drag  coefficient 


2D 


O^ref 

Base  drag  coefficient 
Rotating  band  drag  coefficient 
Viscous  drag  coefficient 
Wave  drag  coefficient 

Skin  friction  factor  for  flat  plate  viscous  flow 
Empirical  constant 

Conversion  factor  between  flat  plate  and  cylindrical  viscous  flow 
Drag  force 

Operational  parameter 


=  k  + 


M, 


0 


Mach  number 
Muzzle  Mach  number 
Operation  parameter 


Ab 

2m 


Radius  of  nose  ogive 

Reynolds  number 

Time  of  flight  of  projectile 

Intercept  of  C^-M  characteristic 

Representative  diameter 
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LIST  OF  SYMBOLS  (Continued) 


Base  diameter 

d^  Rotating  band  diameter 

k  Slope  of  C^-M  characteristic 

Z  Length  of  afterbody  of  projectile 

a 

Length  of  nose  of  projectile 
m  Mass  of  projectile 

v  Velocity  of  projectile 

V  .  Sonic  velocity  at  ambient  conditions 
AV  Velocity  decrement  of  projectile 

p  Ambient  air  density 
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